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ABSTRACT 


~  A  We  consider  the  effect  of  the  differential  rotation  of  the  disk  of 
the  Galaxy  on  magnetic  field  which  penetrates  the  disk.  The  magnetic 
field  will  be  progressively  distorted  from  a  potential  (current-free)  form 
and  will  at  some  stage  become  unstable.  We  expect,  from  knowledge  of 
solar  flares,  that  an  MHD  instability,  a  resistive  instability,  or  a  combin¬ 
ation  of  the  two,  will  result  in  the  release  of  the  excess  magnetic  energy 
and  that  part  of  the  released  energy  will  be  converted  into  heat.  By 
estimating  the  energy  release  and  the  rate  at  which  this  process  will 
occur  and  by  assuming  that  this  energy  input  is  balanced  by  radiation,  we 
obtain  estimates  of  the  parameters  of  the  resulting  plasma.  It  appears 
that  this  process  alone  can  heat  a  galactic  corona  to  temperatures  of  order 

to6  k.  p  OOC  K. 
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1.  INTRODUCTION 


The  diffuse  soft  x-ray  (E  <1  koV)  background  lias  been  attributed  to 
a  coronal-type  plasma  of  T  -  10  K  which  occupies  a  substantial  fraction  of 
the  galactic  volume  (see  I'anaka  &  Bleeker,  1977,  for  a  review).  Models 
for  the  creation  of  this  high-temperature  interstellar  medium  ( ISM)  componen 
through  supernova  explosions  have  been  developed  by  Cox  and  Smith  (1974), 
Smith  (1977),  and  McKee  and  Ostriker  (1977).  Spitzer  (1956),  Shapiro  and 
Field  (1976),  and  Chevalier  and  Oegerle  (1979)  have  analyzed  the  dynamics 
and  energy  balance  of  the  hot  ISM  component;  using  the  supernova  heating 
mechanism,  they  predict  the  formation  of  a  corona  extending  several  kilo- 
parsecs  or  more  above  the  plane  of  the  Galaxy.  We  wish  to  suggest  an 
alternative  process  for  heating  such  a  corona. 

Various  lines  of  evidence  indicate  that  the  Galaxy  is  permeated  by 
a  magnetic  field,  possibly  of  complex  structure.  For  the  same  reason  that, 
if  a  magnetic  field  originates  in  the  solar  photosphere,  it  is  energeti¬ 
cally  favorable  for  the  field  to  extend  far  out  into  the  corona,  we  may 
infer  that,  if  there  is  magnetic  field  in  the  disk  of  the  galaxy,  this 
field  will  also  extend  far  outside  t  lie  disk  into  what  may  be  termed  a 
"coronal"  region. 

In  the  case  of  the  sun,  energy  is  released  Impulsively  at  coronal 
heights  during  flares  (Svestka,  1976;  Sturroek,  1979).  Tt  is  generally 
agreed  that  the  energy  released  during  a  flare  is  stored  above  the  photo¬ 
sphere  in  the  form  of  the  free  energy  of  a  current-carrying  magnetic  field. 
The  influence  of  differential  photospheric  mot  ion  on  magnetic  flux  thread¬ 
ing  the  photosphere  necessarily  leads  to  non-potential  field  configurations 
and  hence  to  free  magnetic  energy. 


Given  the  fact  that  the  disk  of  the  Galaxy  Is  In  a  state  of  differ¬ 
ential  rotation,  and  given  the  reasonable  assumption  that  magnetic  field 
threads  the  disk  of  the  Galaxy  and  extends  out  into  a  galactic  corona, 
the  same  Ingredients  are  present  In  the  Galaxy  as  lead  to  flare  activity 
in  the  case  of  the  sun.  The  purpose  of  this  article  is  to  propose  that 
galactic  flares  provide  a  substantial  energy  Input  to  the  Interstellar 
medium,  and  thus  heat  the  coronal  plasma  responsible  for  the  diffuse  soft 
x-ray  background. 
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II.  ANALYSIS 


This  analysis  will  be  confined  simply  to  order-of-magnitude  estimates 
of  the  effects  of  flare  activity  on  the  galactic  corona.  The  arguments 
and  analysis  will  run  parallel  to  an  earlier  similar  discussion  of  solar 
flares  (Sturrock,  1974).  We  consider  what  is  initially  a  closed  magnetic- 
field  structure  of  mean  strength  B  (gaus.s)  in  the  form  of  a  flux  tube. 
Because  of  the  differential  rotation  of  the  Galaxy,  both  ends  of  such 
a  flux  tube  are  being  twisted  in  the  same  sense,  as  viewed  from  the 
disk  of  the  Galaxy,  but  in  opposite  senses  with  reference  to  the  flux 
tube  itself.  If,  at  some  time,  a  particular  flux  tube  contains  a  potential 
field,  then  some  time  later  the  twisting  of  the  field  will  lead  to  a 
helical  configuration  of  the  magnetic  field.  Unless  the  plasma  stresses 
are  large,  the  field  will  tend  to  adopt  a  force-free  configuration. 

Theoretical  studies  indicate  that  a  force-free  field  becomes  MHD- 
unstable  when  the  distortion  from  the  current-free  form  reaches  a  critical 
level.  Numerical  studies  (Barnes  and  Sturrock,  1972)  indicate  that,  when 
the  distortion  exceeds  the  critical  level,  a  closed  magnetic-field  con¬ 
figuration  is  likely  to  erupt  to  form  an  open  configuration  which  neces- 
arily  involves  one  or  more  current  sheets.  These  sheets  can  then  reconnect 
by  nonlinear  processes  (Petschek,  1964)  arising  from  the  tearing-mode 
instability  (Furth  et  al.,  1963),  returning  the  field  to  its  current-free 
form  and,  in  the  process,  releasing  a  large  fraction  of  the  initial  free 
energy  associated  with  the  coronal  currents. 

Another  somewhat  simpler  possibility  is  that  a  sheared  magnetic-field 
configuration  is  at  some  stage  subject  to  resistive  instabilities  which 
immediately  return  the  field  to  the  current-free  form  (Spicer,  1977). 
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Studies  carried  out  during  the  recent  Skylab  Workshop  on  Solar  Flares  led 
participants  to  the  view  that  both  processes  probably  occur  in  the  sun 
and  are  responsible  for  different  types  of  solar  flares  (Sturrock,  1979). 

It  is  also  possible  (Sturrock  and  Uchida,  1979)  that  the  heating  of  the 
normal  corona  may  be  due  to  the  same  processes  occurring  in  magnetic  fields 
of  such  configuration  that  the  onset  of  Instability  is  non-explosive 
rather  than  explosive  (Sturrock,  1966). 

We  suppose  that,  when  a  galactic  flare  occurs,  energy  in  the  amount 
1  2 

e  ^  B  per  unit  volume  is  released  in  such  a  form  as  to  heat  the  coronal 
plasma.  The  coefficient  e  will  be  the  product  of  two  terms,  and  e^, 
where  is  the  ratio  of  the  magnetic  free  energy  to  the  initial  current- 
free  magnetic  energy,  and  e9  is  the  fraction  of  the  magnetic  free  energy 
which  goes  to  heat  the  plasma.  The  calculations  of  Barnes  and  Sturrock 
(1972)  indicate  that  flare  action  can  occur  after  magnetic  field  lines 
are  twisted  by  half  a  complete  rotation  and  that,  at  that  time,  »  1. 

We  make  the  simple  assumption  that,  during  an  eruption,  one  half  of  the 
free  energy  of  the  force-free  field  survives  as  the  free  energy  associated 
with  the  current  sheet.  We  further  assume  that,  when  reconnection  occurs 
in  the  current  sheet,  one  half  the  energy  goes  into  heat  and  the  rest  into 
mass  motion.  Hence  we  assume  that  e9»l/4  so  that  cm  1/4.  Assuming  that 
the  energy  released  during  one  flare  is  almost  all  radiated  away  by  the 
time  of  the  next  flare,  we  see  that,  directly  after  a  flare, 

3nkT  -  cj^  B2  (2.1) 

where  the  gas  is  assumed  to  be  almost  pure  hydrogen  which  is  almost  fully 
ioni.’.ed  and  n  is  the  density  of  either  electrons  or  protons. 

If  the  plasma  loses  energy  primarily  by  radiation,  we  may  introduce 
the  radiative-cooling  time  scale  T  defined  by 

K 
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n 


(2.2) 


-  3nkT  t"1 

K 

-3  -1 

where  n  (erg  cm  s  )  is  the  radiative  rate.  We  adopt  for  n  the  form 

n  -  Tn2TY  (2.3) 

given  by  McKee  and  Cowle  (1977)  in  the  range  10^ < T <  10^ ‘ 6 ,  where  the 

—  18  2 

radiation  rate  Is  well  represented  by  the  above  equation  with  T  -  10 
and  Y  -  -3/5. 

The  effect  of  differential  rotation  of  the  Galaxy  is  such  that  any 
element  of  area  exhibits  an  "orbital"  rotation  and,  in  addition,  a  "planet¬ 
ary"  rotation  which  may  be  characterized  by  a  vorticity  given  by 


to  = 


(2.4) 


where  fl(r)  is  the  angular  velocity  of  the  Galaxy  at  radius  r.  Hence  a 

small  flux  tube  which  emerges  from  one  location  in  the  disk  and  returns 

to  another  location  at  about  the  same  radius  will  be  twisted  at  both  ends 

in  such  a  way  that  the  flux  tube  is  given  a  helical  twist.  We  expect  that 

the  tube  will  become  MHD  unstable  when  the  twist  reaches  a  critical  value. 

Numerical  studies  (Barnes  and  Sturrock,  1972)  indicate  that  this  value 

corresponds  to  half  a  complete  rotation  which  would  result  from  a  rotation 

of  n/2  at  each  end.  Hence  we  may  estimate  the  time  between  flares  t  (s) 

F 

from 


to  t  s  n/2 .  (2.5) 

r 

We  find  from  Mihalas  and  Routly  (1968)  that,  over  the  range  of  radii 

3  kpc  to  13  kpc,  to  has  a  value  close  to  30  km  s  ^kpc  ^  or  10 

15  2 

Hence  we  obtain  the  estimate  r  -  10 

F 

If  the  corona  plasma  reaches  an  approximate  steady  state  in  which 
heating  by  flares  is  balanced  by  cooling  by  radiation,  then 
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(2.6) 


XR=  TF  =  T  • 


We  obtain  from  (2.1),  (2.2),  and  (2.3),  the  following  estimates  of  the 
time-average  values  of  the  density  and  temperature: 


»  -  109-3  B1W13  t-5/13 


T  -  ID4'1  B10/13  T5'13  , 


(2.7) 


which  for  the  above  estimate  of  T  becomes 

n-  103-5B16/13  , 

T  -  IQ9’9  B10/13  . 


(2.8) 


If  the  magnetic  field,  which  contains  the  coronal  plasma,  extends  out  to 
a  distance  L  (cm)  from  the  disk,  then  the  emission  measure  H,  defined  by 

H  -  n2  L,  (2.9) 


will  have  the  value 

H  =  IQ7'0  B32/13  L  . 


(2.10) 
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III.  DISCl'SSION 


Analysis  of  the  soft  x-ray  background  Indicates  that  the  plasma  respon¬ 
sible  for  this  emission  probably  has  a  complex  thermal  structure  (see, 
e.g.  Cash  et  al.,  1976;  Burstein  et  al.,  1977;  Stern  and  Bowver,  1979). 

Nevertheless,  for  the  order-of-magnitude  analysis  presented  here,  we  may 

5  8 

adopt  as  typical  coronal  parameters  a  temperature  of  ~  10  K  and  an  emis¬ 
sion  measure  of  10^“  cm  ^  (5.10  ^  pc  cm  We  see  from  equation  (2.8) 

that  this  temperature  is  consistent  with  the  present  model  if  the  mean  magnetic 

-5  3 

field  strength  in  the  local  region  of  the  galactic  corona  is  -  10  ’  gauss. 

For  comparison,  we  note  that  Badhwarand  Stephens  (1977),  on  the  basis  of 
an  analysis  of  the  hydrostatic  equilibrium  of  the  gas-field  system  in  the 
solar  neighborhood,  conclude  that  the  magnetic  field  strength  in  the 
neighborhood  of  the  sun  is  about  (5-6)  x  10  h  gauss,  and  that  Ruzmajkin 
and  Sokolov  (1977),  on  the  basis  of  an  analysis  of  pulsar  data,  conclude 
that  the  local  magnetic  field  strength  is  (2.1  *  1,1)  x  10  ^  gauss.  Our 
'heoretical  estimate  of  the  required  magnetic  field  strength  agrees  with 
these  two  (observat lonal lv  derived!  estimates  as  well  as  they  agree  with 
each  other. 

We  now  see  also  from  equation  (2.8)  that,  with  the  above  estimate  of 

the  magnetic  field  strength,  we  expect  the  electron  density  in  the  galactic 

-3.0  -3 

corona  to  be  about  10  ‘  cm  .  We  also  see  from  equation  (2.9)  that  the 

16  2 

inferred  emission  measure  H  = 10  '“is  obtained  if  the  corona  extends  out 

22  2 

to  a  distance  L « 10“  '  ,  i.e.  5kpc. 

It  is  necessary  to  compare  the  above  estimate  of  L  with  the  scale 
height  h  (cm)  of  plasma  in  the  galactic  corona.  We  use 


if  the  plasma  is 


and  note  that  the  mean  mass  m  (g)  has  the  value  10 
mainly  fully  ionized  hydrogen,  and  that  the  gravitational  acceleration 
g  (cm  s  is  approximately  10  '  in  the  range  1  kpc  to  5  kpc  of  the 

galactic  plane  (All^n,  1973;  p.  250)  ,  so  that,  numerically, 

h  » 1016,3  T.  (3.2) 

5  8  22  1 

We  find  that,  for  the  estimated  temperature  T«10"  ,  h  *»  10  .  Hence, 

if  the  value  of  L  is  in  fact  set  by  the  extent  of  the  magnetic  flux  tubes 
which  contain  the  plasma  in  the  local  region  of  the  galactic  corona,  the 
scale  height  is  sufficiently  large  to  fill  the  flux  tubes  with  plasma. 

If,  on  the  other  hand,  the  magnetic  flux  tubes  extend  well  beyond  5  kpc, 
then  the  extent  of  the  coronal  plasma  may  be  attributed  to  the  scale  height. 

If  the  latter  is  the  correct  interpretation,  we  may  use  (3.2)  in  (2.10) 
to  obtain  an  expression  for  H  in  terms  of  B  only.  Then  (as  we  have  seen) 
our  choice  of  B  fits  both  T  and  H.  However,  our  initial  equations  [(2.1) 
and  (2.2)] are  based  on  the  assumption  that  n  and  T  are  uniform  for  any 
flux  tube,  and  this  assumption  must  be  relaxed  if  the  scale  height  is  less 
than  the  height  of  the  flux  tubes. 

The  fact  that  the  present  model  gives  approximate  order-or-magnitude 

agreement  with  observational  data  does  not  establish  its  validity;  it  simply 

indicates  that  more  detailed  investigation  should  be  pursued.  In  particular, 

it  is  necessary  to  relax  the  "steady-state"  assumption  and  to  investigate 

the  time  development  of  heating,  cooling  and  mass  flow.  The  condition 

(2.1)  is  valid  directly  after  a  sudden  release  of  energy,  whereas  condition 

(2.6)  is  valid  for  a  long-term  average  of  the  coronal  conditions.  If 

were  in  fact  short  compared  with  T  , then  energy  released  during  a  flare 

r 

would  rapidly  radiate  away  so  that  the  time-average  temperature  would  be 
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less  than  the  Immediate  post-flare  temperature.  If,  on  the  other  hand, 

t  were  large  compared  to  T„,  then  energy  released  during  flares  would  accumu- 
K  r 

late  with  the  possible  result  that  the  plasma  stress  eventually  exceeds 

the  possible  magnetic  stress.  This  would  result  in  an  opening  of  the 

magnetic  field  with  loss  of  the  coronal  plasma. 

If  T  (s)  is  the  time  scale  for  the  release  of  energy  in  a  flux  tube, 

E 

then  (Sturrock,  1974) 


(3.3) 


where  R  (cm)  is  the  radius  of  the  flux  tube  under  consideration  and 

v  (cm  s  )  is  the  reconnection  rate.  Following  Petschek  (1964),  we  esti- 
K 

mate  that 


v  =  10-1  v  (3.4) 

R  A 

where  v  is  the  Alfven  speed.  Since  the  mass  density  of  fully  ionize^ 

A 

hydrogen  is  2nm,  where  m  is  the  mean  particle  mass, 

v  =  B(8Tinm)~1/2  *  1011,3  B  n"1/2  (3.5) 

A 


-5  3 

If  we  consider  the  value  B  -  10  ",  we  find  that,  for  R  in  the  range 

10“A'3  to  10““'“  (1  kpc  to  5  kpc),  t  is  in  the  range  103^’^  to  lO33'^. 

h 

It  appears,  therefore,  that  Tg  is  comparable  with  Ty  This  suggests  that 
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the  corona  will  be  approximately  in  a  steady  state  so  that  the  heating  is 
more  like  the  steady  heating  of  an  active  region  rather  than  the  impulsive 
heating  of  a  flare. 


We  should  also  examine  the  rate  of  energy  loss  through  thermal  conduc¬ 
tion  to  the  galactic  disk,  to  compare  with  the  energy  loss  rate  by  radia¬ 
tion.  If  T  (s)  is  the  time  scale  for  cooling  of  the  galactic  plasma  by 
H 

thermal  conduction,  then  x  may  be  estimated  from 

n 

3  nkTL  =  F  xu,  (3.8) 

H 

-2  -1 

where  F  (erg  cm  s  )  is  the  heat  flux  along  the  magnetic  field.  This 
may  be  estimated  from 

F  =  K  T  L-1=  10" 6  T?/2  L_1,  (3.9) 

if  we  use  the  Spitzer  (1962)  form  of  the  thermal  conduction  coefficient. 
Hence  we  obtain  the  estimate 


l2„t-5/2. 


<3.10) 


For  the  values  L  =  1022'2,  n  =  10  ^  and  T  =  10^’^,  we  obtain  x^  =  102^’^, 
giving  a  cooling  time  of  1010y.  This  is  much  longer  than  the  radiation 

cooling  time,  confirming  that  radiation  is  the  dominant  mechanism  for  energy 
loss.  On  the  other  hand,  X  is  comparable  with  the  age  of  the  Galaxy,  so 
that  the  energy  flow  from  the  corona  to  the  galactic  disk  is  sufficient  to 
produce  the  "evaporation"  (Neupert,  1968;  Antiochos  and  Sturrock,  1978) 


necessary  to  provide  the  coronal  plasma. 


It  has  been  known  for  some  time  that  there  are  high-velocity  clouds 
(HVC)  and  intermediate-velocity  clouds  (IVC)  of  neutral  hydrogen  outside 
the  galactic  plane.  Verschuur  (1975),  in  reviewing  this  topic,  suggests 
that  there  is  more  than  one  type  of  cloud  and  notes  several  possible  inter- 


pretatlons.  Shapiro  and  Field  (1076)  have  developed  a  "galactic  fountain" 
model  involving  an  upward  streaming  hot  gas  which  cools  and  condenses  into 
a  cool  component.  Our  comparison  between  the  galactic  corona  and  the  post  - 
flare  plasma  in  the  sun  suggests  that  some  of  these  clouds  may  he  the 
galactic  analogues  of  "coronal  rain"  which  often  occurs  In  the  loop  promi¬ 
nence  systems  formed  by  some  large  solar  flares  (Tandberg-Hanssen ,  19741 . 

It  has  been  shown  by  Goldsmith  (1*1711  and  by  Ant lochos  (1976)  that  the 
latter  may  he  understood  as  the  result  of  thermal  instability  in  the  hot 
dense  plasma  which  forms  in  the  coronal  volume  of  an  active  region  as  the 
result  of  a  flare. 

Following  Field  (1965)  and  Ant lochos  (1*176),  we  note  that  the  condi¬ 
tion  for  thermal  Instability  is  that  the  wave  number  k  (cm  *)  of  the  small- 
amplitude  wave  should  satisfy  the  inequality 


k‘ 


j 

k‘ 

c 


an 

On 


(3.11) 


On  using  the  form  (2.3)  for  n  and  the  Spitzer  form  for  the  thermal  con¬ 
ductivity  k,  we  find  that 


3  1  /  •>  l/->  -v  /  2  —  7/4 

k  -  10  (2  - 1)  '  r  nT’  , 

c 


(3.12) 


With  the  values  of  >  and  I'  appropriate  to  the  range  10^  T  <  I0^*h,  this 


becomes 


k 


10 


■5.0 


-41/20 


(3.13) 


—  3  *>  8  —1  20  8 

For  the  values  n  -  10  ,  T  -  10  *  ,  we  find  that  k^  »  10  *  .  Since 

this  estimate  of  k'  Is  considerably  less  than  our  estimate  of  1.,  we  expect 

that  thermal  Instability  will  lead  to  the  formation  of  condensations 

with  dimensions  (along  the  magnetic  field)  of  200  pc  or  more.  When  formed. 
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